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Abstract 



We formulate theory of interacting scalar field on the fuzzy sphere as a random matrix 
model. We then analyze the expectation values of observables of the theory in the large 
N limit and we demonstrate that the eigenvalue distribution of the matrix M remains the 
polynomially deformed Wigner semicircle. We also compute distributions involving the 
matrix Laplacian of M and we show that the correlation between the eigenvalues of these 
two is different from the free field case. 



1 Introduction 



Matrix models have established a firm place in the modern physics. Starting with pioneering 
work of Wigner in description of the spectrum of heavy nuclei, they have emerged in areas of 
string theory as finite approximations to Riemann surfaces [2], condensed matter system as 
the Calogero model [3] or chaos in quantum systems [Ij. Fuzzy spaces are non-commutative 
spaces with a finite dimensional underlying Hilbert space |5]. It is possible to describe 
them by finite dimensional matrices and they become their commutative counterpart in the 
limit of these matrices being very large. The physical motivation of using fuzzy spaces is 
in regulating the divergences without breaking the isometrics of the space-time, which is 
especially appealing to quantum gravity considerations [6] . Fuzzy spaces also arise as brane 
solutions in string theory and in M-theory |7j. 

In [8j, it was shown that matrix models and fuzzy spaces are related. One can generalize 
the standard Gaussian matrix ensembles by addition of matrix Laplacian. This procedure 
is motivated by the the scalar field theories on fuzzy spaces and introduces new observables 
involving matrix derivatives. Averages of some of these were considered. This was then used 
to show, that the eigenvalue distributions of the matrix and its Laplacian are correlated and 
the joint distribution was computed. From physical point of view, such ensembles provide 
new computational tools for the fuzzy theories, for example to compute distributions of 
various observables. 

In this paper, we present some further work along these lines. In the first part, we 
compute expectation values and derive the distributions of more general observables in 
such defined random matrix ensembles. This shows that this framework is very robust 
and can be used beyond simple observables of the form M^B^, where B is the matrix 
Laplacian of M. We also show, that there is a different scaling of the terms in the action, 
under which the contribution of the mass term survives the limit of very large matrices. 
Moreover, the correlation between the distributions of M and B then remains finite for any 
form of the kinetic term in this limit. In the second part, we concentrate on the matrix 
ensemble motivated by an interacting field theory. We show that, as was the case before, 
the distribution of eigenvalues of M is altered only by a rescaling of the variable. We then 
investigate the properties of observables involving B. We show that the distribution of the 
eigenvalues of B and the joint distribution of M and B are connected and we compute each 
fo these up to the second order in the MB correlation. We also show that the interaction 
brings some new features into the problem, which were not present in the free case. 

Some aspects of this problem were approached by other authors from different points of 
view. Authors of [9] treated the Laplacian term as a perturbation and integrated out the 
angular degrees of freedom. The eigenvalue problem was then solved using the standard 
methods. After the corresponding approximation is made, our results presented here agree 
with the results computed using this method in |1UJ. In the kinetic term and added 
interaction term are treated exactly, but the eigenvalue distribution is said to have the 
same destiny as in the free case and different aspects of the results are considered. Here, 
we present more complete calculation of the distribution and we also give calculation of 
expectation values of observables involving the Laplacian. 



2 



2 Results for the free theory 



We will consider the Euclidean theory of a real scalar field on the fuzzy sphere governed by 
the general kinetic term action 

5o(M) = ^Tr (M/CM) + ^fi^Tr (M^) . (1) 

We will introduce interaction terms later and for now we will consider only the free theory. 
The standard Laplacian kinetic term is given by /CM = [L^, [Lq,,M]], where Lq, are the 
generators of the N dimensional representation of SU (2). We denote /CM = B. We further 
introduce a basis for fields on the fuzzy sphere in terms of N x N matrices 

, / = 0,l,...,iV-l , m = -l,-l + l,... ,1-1,1 , (2) 

normalized as Tr (T^ T^, ) = (5'' 6m-m' ■ We can expand the matrix M in terms of this basis 
as 

M = Y,ciTl (3) 

l,m 

The non-interacting action is diagonal in this basis and the correlator of two components of 
M is (cjj^ ^m'/ — ^ ^mm'GiV) where the propagator G{V) depends on the form of the kinetic 
term. For example for the standard kinetic term it is G{1) = l/(^^ + l{l + 1)). 

We define the following two-point functions 

{{MM)^j) = f5^J , {{BB),j) = g5^J , {{MB)^j) = h6^j. (4) 

Also / = (Tr (MM)) /N and similarly for g and h. In fact, these three correlators are 
all the information we need about the theory and we do not need to know the form of /C 
explicitly. This also means that with a proprer choice of /C, we could work with the theory 
on different fuzzy spaces. 



2.1 Previous results 

Here, we very briefly summarize the results of [8], which we are going to use later. In this 
section, we will normalize the distributions of M and B to have unit radius. This eliminates 
some cumbersome factors of 2 in final formulas. In the sections to follow, we will however 
change this normalization to eliminate factors of two from equations we will work with and 
to follow the standard convention. This should be kept in mind when comparing results 
from these sections [2] and HI 

We define the normalized correlators 



m,b 




M 



B 



(5) 
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We can write down, in the large N limit, the following recursion rules for them 

m-2 6-1 
^Wm,b = Yl ^P,oWm-2-p,b + lY^ Wm-l,b-l-pWo ,p m > 1, 

p=0 p=0 

6-2 m-1 
4W-„,6 = J2^^'P^rn,b-p-2 + l Y.^"^-^-P^b-lWp,o b > 1, (6) 

p=0 p=0 



where 7 = h/^/fg. These are results of explicit Wick contractions of matrices and the 
planarity of corresponding diagrams. Defining the generating function 

00 

<p{t,s)= J2 Wm,bt"'s' (7) 

m,6=0 

these recursion rules become equations for (j){t, s), which can be solved as 

0(i,O) = <^(t) = 2i^^^^H? = _|_ , 0(O,s) = </.(s) , (8) 

Inverting this generating function we arrive at the final formula for the distribution function 

1-72 

p{x, y) = p{x)p{y) 2\2 TTT, — 2^ , a 2f 2 ±_ 2^ ' (1^) 

(1 — 7"=)^ — 47(1 + 7^ja;y + 47^(x^ + y^) 

where p{x) = 2\/l — z^/vr is the Wigner semicircle distribution generated by (j){t). This 
distribution is positive for 7^ < 1 and in the limit 7 — > ±1 becomes p{x)6{x — y), so the 
two matrices are completely correlated/anti-correlated, as expected. 

Therefore the distribution of eigenvalues of the unnormalized M has radius 2^/J, distri- 
bution of eigenvalues of B has radius and the correlation between the two is given by 
the above formula. 



2.2 Joint distribution of eigenvalues of three matrices 

Since the two point functions of the matrices M and B were the only relevant quantities for 
computation of the joint distribution of the previous section, if we define 'linear' matrices 
A as 

Ai,A2 = 5;4/i,2(/,m)r^, (11) 

m,l 

the whole procedure will go through basically intact and we recover the same results. 
Namely the joint distribution of the eigenvalues of Ai and A2 is going to be given by 

1 - 7^ 

p{x, y) = p{x)p{y) 2^^^ — A 2 f 2^ 2^ (^^) 

(1 - 112) - 4712(1 + 7t2)xy - 47f2(a;^ + y^) 



4 



and 712 = {Tr{AiA2)) /y^O]^, where aj = (TrAfy The distributions of eigenvalues of 
matrices Ai, A2 have radius 2^ai 2- In terms of the original definition we have 



N-l I 



(13) 



1=0 



m=~l 



We now turn to computation of the joint distribution of three matrices of the form (11), 
i.e. we look for p{x, y, z) such that 



A-i 



A2 Yf Ai 



dxdydz x"'y^z'^p{x,y, z). (14) 



Our approach will be again to write down recursion rules for the moments Wa,b,c coming 
from the Wick contractions. There are three different types of contractions 



Tr ( AIA^A^ 



Tr(Ar2))aiiV(Tr(^«^^ 



\b AC 



(Tr (^r^^^)) y^7i2iV (Tr Ut^-M^ 



Tr 1 V^7i3iV (Tr (A^)) . 



(15) 



Each contraction splits the diagram into two parts and planarity condition forbids contract- 
ing matrices in different parts. Summing over all possible contractions we obtain 



a- 2 



6-1 



4Wa,6,c = E ^P,0,0Wa-p-2,b,c + 7l2 ^ Wa-l,pflWo^b-p-l,c + 



c-1 



+ 713 — l,b,c—p—l 

p=0 



a> 1. 



(16) 



This rule holds only for a > 1, because for a = there is no matrix Ai to do the contraction 
with. The quantities Wo,fe,c are moments of two point distribution of matrices A2 and ^3 
and are given hy In this recursion rule they enter as the initial condition. 

Considering contractions of the other groups of matrices we obtain two more equations 

b-2 a-l 



b,c 



Y.W^,p,^Wa,b -p-2,c + 712 -p-l,6-l,c + 



p=0 



p=0 



c-1 



+ 723 Y.Wafi,pWo,b -l,c-p-l 6 > 1, 



(17) 



p=0 



a-l 



Wofi,pWa,b,c-p-2 + 713 E ^0,p,0l^a-p-l,fe,c-l + 



b-1 



+ 723 — 1— p,c— 1 

p=0 



c > 1. 



(18) 
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We solve these by defining the generating function 

s, = (19) 

a,b,c 

and by rewriting the recursion relations as equations for (f)(t,s,u) 

4(^(/>(t, s, u) - (/>(0, s, n)) = 0, 0)(p{t, s, u) + 7i2ts</'(t, s, 0)0(0, s, u) + 

+ 7i3tu0(O,O,n)</.(t,s,'u), (20) 

A(j){t, s, u) - <j){t, 0, n)) = sV(0, s, O)0(t, s, u) + 7i2ts</'(i, 0, 0)(?;)(t, s, n) + 

+ 723S^^0(i,O,n)0(O,t,n), (21) 

4^(/)(t, s, -u) - (j){0, s, u)^ = u'^(j){0, 0, u)4>{t, s, u) + 7i3tii(/)(0, s, 0)(/)(t, s, u) + 

+ 723Stx0(i,s,O)</.(O,t,n). (22) 

We require 0(0, 0, 0) = Wo,o,o = 1 as a normalization condition. Now, choosing appropriate 
variables to be zero, we can solve these equations. For example in the first relation, setting 
s = li = we get 

40(t,O,O) -4 = t202(t,O,O), (23) 
which is the same equation we have arrived at in the two matrix case and has solution 

0(t,O,O) = 0(t) = ^^--^==. (24) 
Similarly for the case of 0(0, s, 0) and 0(0, 0, u). 



Setting n = in (20) we get equation for 0(t,s,O) which solves again for the formula 
obtained in the previous section 

4 - t^0(t) - 7ts0(s) ^ - 7ts0(s) 1 - |7i2ts0(t)0(s) 



We could obtain the same quantity from (21 ) by setting u = and this obviously yields the 
same result. Next, we set s = in the first equation and obtain 

0(^,0, n)= ^ ^ (26) 

1 - -ili?,tu<t>(t)(j){u) 

and the same way for the rest of the functions. 

Plugging these into one of the original equations we get the final formula for the gener- 
ating function 

/ , N 0(g)0(t)0(n) 

1 - ist7i20(s)0(t) 1 - ^tu7i30(t)0(n) 1 - isu7230(s)0(n) 



With no surprise, this is the formula we obtain from either of the three equations. 
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We will now describe a general method that can be used to invert generating functions 
of this form to obtain the corresponding distribution. We present the proof in the appendix 
[A} If the generating function is expressed as 

(^(tl, ...,tn) = fih^ih), tn(l){tn))Hh) ■ ■ ■ 0(tn), (28) 

then the corresponding distribution is given by 

p(xi , . . . , x„) = ^ £i . . . enF{e''^'^ e*^"^" ) , (29) 

£i=±l 

where 



F{zi, • • • , ^n) = I n ^ 1 • • • ' (30) 



and Xji — COS0J7,. Using this for the generating function »2,7y gives a very complicated 
formula of the form 



P{x, y, z) = p{x)p{y)p{z) X p3(x, y, z). (31) 

Explicit formula for p3{x,y, z) is given in the appendix [B| 

However couple of important observations can be made. In the case 712,23,13 = the 
factor becomes 1, in the case of two of the three 7's vanishing the factor becomes the 
appropriate function to give p{x,y,z) = p{x,y) x p{z), if 712,23,31 — )• 1 we get p{x,y,z) — )• 
p{x)5{x - y)5{x - z) and if 712,23 -1,723 1 we getp(x,y,z) p{x)5{x + y)5{x + z), 
i.e. fully correlated or anti-correlated distributions as expected and finally in the case of 
712 = 713 = 7,723 1 we obtain p{x,y,z) p{x,y)5{y - z). 



2.3 MBMB joint distribution 

To get the four point joint distribution of matrices MBMB, we need to compute the 
following quantities 

Wa,h,c,d = ^ ^Tr 

Using the same explicit Wick contractions and planarity of the diagrams we find the fol- 
lowing large N recursion rule 

a-2 6-1 

,0,0,0 W^a-p-2,6,c,d + 7^] ^0,p,0,0^a-l,6-p-l,c,d + 
p=0 p=0 
c-1 6-1 

-l,0,c-p-l,d + 7 ^ ^a-l,0,0,p^0,6,c,d-p-l (33) 
p=0 p=0 



M Y 

V7y 



B 



VV7 



B 



(32) 
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for a > 1. We have considered all the possible contractions of the first matrix of the 
part. Introducing the generating function 

cl){t,s,u,v)= Wa,b,c,dt''sW (34) 

a,b,c,d 

this becomes 

i(^4>it, s, u, v) - 0(0, s, u, v)^ = t'^(f){t, 0, 0, O)0(t, s, u, v) + 7ts(/)(0, s, 0, 0)0(t, s, u, v) + 

+ tu(/>(0, s, u, O)0(t, 0, u, v) + 'ytv(j){t, 0, 0, u)(/>(0, s, u, f X-35) 

Now if we consider contractions of matrices from other parts of the diagram, we get three 
more equations for the generating function, namely 

4^(/)(t, s, n, v) — <t){t, 0, u, v)^ = s^0(O, s, 0, 0)</>(t, s, n, + jsu(j){0, 0, n, O)0(t, s, u, v) + 

+ si;0(0, 0, u, v)0(i, s, 0, v) + 'yst(j){t, s, 0, 0)(p{t, 0, u, u), (36) 

4^(/)(t, s, u, v) — (j){t, s, 0, v)^ = u'^(j){0, 0, u, 0)(j){t, s, u, v) + 'yuv(p{0, 0, 0, v)<t){t, s, u, v) + 

+ ut<l){t, 0, 0, v)(j){t, s, u, 0) + 7us</>(0, s, n, 0)(/)(t, s, 0, u ),(37) 

4(^0(t, s, M, w) - (j){t, s, M, 0)^ = 1-^0(0, 0, 0, ■u)(/'(t, s, ti, w) + ^vt<j){t, 0, 0, 0)(/)(t, s, u, v) + 

+ vs(j){t, s, 0, 0)0(0, s, u, u) + jvu(j){0, 0, M, v)(l){t, s, u, 0)(38) 

We again have 0(0,0,0,0) = VFo,o,o,o = 1- Then, choosing appropriate variables to be zero, 
we can solve these for the generating function pretty much the same way we did in the case 
of three matrices, with the same results for (p{t, 0, 0, 0), 0(t, s, 0, 0), 0(i, s, u, 0) and the rest 
of the combinations. We can make the notation more compact by defining 

i;{t,s) = l-^jtsts^{t)(P{s), (39) 

where 'jts is the correlation parameter between the matrices corresponding to the variables 
t and s. This way 

0(r)0((j) 

(P[T,cr} = — -, (40 

.f N 0(r)0(o-)0(p) 

where r, o", p are any of t, s, u, v. Finally (32) becomes 

0(t)0(s)0(u)0(^;) 



(j){t,S,U,v) 



ip{t, s)0(t, v)4'{u, v)ip{s, u)ip{s, v)^l){t, u) 



lb J 

Using the method described in detail in the previous section, this yields the four point 
distribution of the form 

p{x, y, z, w) = p{x)p{y)p{z)p{w) x pi{x, y, z, w), (43) 

with explicit formula for p4{x, y, z, w) is given in the appendix [B| Here we just observe that 
the factor is 1 when 7 = 0, in the case of 7 — )• 1 we get p{x,y, z,w) — )■ p{x)6{x — y)S{z — 
w)5{x — w) and in the case of 7 — )• — 1 p{x, y, z, w) — )■ p{x)5{x + y)S{z + w)5{x + w), i.e. the 
variables correctly correlate or anti-correlate. 
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3 Mass rescaling and the correlation 7 



The large N structure of correlation parameter 7 was discussed in [8]. It was shown that if 
G{1) goes like I" for large /, correlation 7 can vanish or tend to a constant value in the large 
limit, depending on the value a. Especially for the Laplacian kinetic term, i.e. a = —2 
the correlation vanishes as 1/ log A^. In this section, we show that there is a different scaling 
of the kinetic and the mass terms of the action and that upon a rescaling of /x correlation 
7 is always finite. 



The MM correlator is given by 



N-l 



/= l(Tr (MM)) = lj; (2/ + 1)G(/). 



(44) 



In the case of Laplacian kinetic term, the large limit of this expression is 



/ 



dx 



2Nx + l 



Nx{Nx + 1) + 



(45) 



which has the advertised log N/N behavior which leads to logarithmic vanishing of 7. If we 
however rescale the mass ^u^ — t- N'^fl'^, we will find that / now depends polynomially on N. 
Namely 



/ 



dx- 



2Nx + l 



/o Nx{Nx + 1) + iV2^2 
Similar calculation then yields also 



1 , ^ 1 



(46) 



and finally finite 



h 
,9 



N 



N-" 



1 - /iMog ( 1 + 



^-A^ + /i^log(l + ^ 



(47) 
(48) 



l-/i2log(l + 



logfl + J. 



(49) 



Note that this is finite in the limit of very large fl and tends to v3/2, which is the correlation 
in the case of no kinetic term. 

The same line of attack works also in the case of a general kinetic term /C. We have 
already assumed, that /CT^ depends only on / and lets assume that for large I, this is 
proportional to 1°. We therefore need to rescale fj, — >• N^fs^, so that for large /, the mass 
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term in the propagator does not get suppressed. The same procedure as before then yields 



/ 
h 



rl—a 



N 
N 



1 



1-2^1 

2 



2 2 

a a 
2 2 

-,i;i + -;- 

a a 



/i2 



2 + a 



2 2 

a a 



(50) 
(51) 
(52) 



where 2F1 is the ordinary hypergeometric function. These clearly give a finite 7. 



Therefore an appropriate rescaling of the mass makes the mass term always unsuppressed 
and keeps the correlation parameter 7 finite for any kinetic term. This is important when 
one sets to apply these results in the field theory, since then we want all the terms to survive 
the large limit. 



4 Interacting theory 

We are now ready to introduce the interaction the free action ([T]). We will consider a quartic 
interaction potential 

5i„t = ^TV(M4) , ~g = g/N. (53) 



As mentioned in the introduction, the case without the kinetic term is well known |12j.|13j 
and the result is a polynomial correction to the Wigner semicircle distribution. 

Expanding the interaction part in power series in g yields for an average of an observable 
0(M) 

(0) = ^ ^ / (iMe-^°(*^)0(M)TV {M^f , (54) 



a=0 



where 



Z = f2 / dMe-^«(*^)Tr (M^)" . (55) 

So we see that evaluating the average in the interacting theory can be done using averages 
of the free theory. We just need to pick correct diagrams that contribute to the expectation 
value on the RHS of ( |54[ ). The diagrams containing vacuum bubbles, i.e. parts, where some 
of the vertexes from Tr (M^)'* contract only among themselves, will be canceled by the 1/Z 
factor, therefore we can write 



_-~9r 

a=Q 

where the subscript will indicate that we consider only diagrams that do not contain dis- 
connected vacuum bubbles and that the contractions are to be taken using the free theory 
measure. 
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Figure 1: Two types of terms contributing to (59) 



Now we also see the motivation for the 1 /N factor in the definition of the coupling con- 
stant g. Since each trace in (|/(M)Tr (M*^)"^ raises large dependence of this expression 
by one, expectation values in the previous sum are going to be all of the same order and all 
the terms will contribute in the large N limit. 

Before we proceed with computation of the eigenvalue distribution of M, let us stress one 



point. The contractions in (56) are done using the free measure. But we have already seen 
that in the free theory, the kinetic term only rescaled the radius of the original distribution. 
And therefore we expect the same in the interacting case, namely that all the distributions 
of the M^-theory with no kinetic term will survive also in the full theory, only with a 
rescaled variable. 



4.1 Eigenvalue distribution of the matrix M 

To compute the distribution of eigenvalues of M we need to compute the moment generating 
function 4>i{t). The potential is even and therefore odd moments vanish and we need to 
compute (Tr M^™). From (56) we see that we need to investigate quantities of the form 



2m 



Expressions 



1 



Tr 





Tr 








Wf) 







a times 



Tr 



M 

7j 



(57) 



O,con 



2m 



a=0 



-F, 



2m 



(58) 



are then going to be finite and will give 2m point correlators of M's in the interacting 
theory. Let us stress again that in this section, recursion rules for the free correlators 
change a little, due to a different normalization of the distribution of M. We will now 
write down the recursion rules for F^^ in a very similar fashion we did in previous cases. 



11 



The contributing planar diagrams are going to be of the form of 2m points on a circle 
with four-point vertexes inside this circle. The number of vertexes is a and the point and 
vertexes connect with no lines intersecting. If we look at one point on the circle, there are 
two different types of contractions we can make. With a vertex, three legs of the vertex 
become effectively three new points on the circle, which now has 2m + 2 points, but there 
are only a — 1 vertexes left and there are 4a different legs that we can connect the first 
point to. Or with a different point, the circle splits into two circles, each new circle can 
have different number of vertexes in it, one circle a and the other a — a, with a proper 



combinatorial factor. This procedure is illustrated in the figure 4.1 We need to be a little 
careful. If there are no points on any circle to connect with, we can not have any vertexes 
in this circle, as this would produce a non-connected diagram, therefore we need to set 



F° = for any a 7^ and = 1. 



Summing over all possible contractions, we obtain 



(Tr (m2-) [Tr (M^)]") = 4a/iv(Tr (m2-+2) [Tr (M^)]""') + 

ma y K 

( J (Tr (m2(p-i)) [Tr(M4)]") (tt (^M^^^'-P^) [TV(M^)]""") (59) 

P=l Q = ^ ^ 



and using the definition (57) this becomes 



Tp2{p-1) p2(m-p) 

a\{a-ay. " ' 



(60) 



p=l a=0 

Immediately we note, that if we define F^"^ = F^"^/a\ this simplifies into 

m—l a 

p2m ^ -p2{m+l) , p2p p2(m-l-p) 

p=0 a=0 

Multiply the recurrence relation by {—g)"" and sum over a = 1 to oo to obtain 



(61) 



m— 1 



F, 



-'2m 



2m — ^0 



(m+l) 

+ 2^ ^F2pF2 

(m-l-p) -?^0 -^0 

p=0 



p2pp2(m-l-p) 



(62) 



We know from the free theory that Fq = Cn and using the identity for the Catalan numbers 



Cn = ^p=Q CpCn-i-p the extra terms cancel. Note that this is consequence of the free case 
recursion rules, and similar terms will cancel for this reason also later. We therefore get 



^2(m+l) 



4<7 



m— 1 



2p-C'2(m-l-p) 



2m 



p=0 



(63) 



This expression holds however only for m > 1, so we need to specify Fq and F2. From the 
definition of F^ it is clear that Fq = 1, consistent with the normalization of the distribution. 
F2 is essentially the dressed propagator of the interacting theory and for a while lets go 
further without specifying it. 
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We now define the moment generating function 4>i (t) = ^ t"^Fm = t^™'P2m , multiply 
the previous formula by ^^(m+i) sum over m = 1 to oo to obtain 

_i_t'F2 = j^ [tUiit) - t\Mt) - 1)] (64) 



or 



which gives 



tUt - (45 + n ct>i + (4(7 + + 4<7t'F2) = (65) 



nit) = . (66) 



So we are left to specify the two point function F2. 

At this point, we are going to take F2 to be the expression obtained by the standard 
methods |12j.|13j. This formula is explicitly 

where a? = (\/l + i&g — l)/2Ag. This might seem that we are assuming something we are 
trying to prove, but this is not the case. We assume that the initial condition we are about 



to use in (63) is the way we expect it to be, i.e. the same as for the case of no kinetic term. 
We then compute the generating function and if it turns out to be the same, we conclude 
that this assumption leads to all the correlators being of the form as in the theory with no 
kinetic term. 



Before we show that this is indeed the case, let us go back to (60) for a while. We can 
use this recursion rule to generate a lot of terms -F^™ and after some trial and error using 
the integer factorization of the terms we can guess the formula for in the following form 

p.2m ^ . 2« ( ^ ( {m + 2a-l)l \ 

\m\{m - 1)1 J \{m + a + l)\a\ J ' ^ ' 

And not surprisingly this is indeed the formula of g expansion of the 2m--point correlator 
given in the appendix of [12]. It would be interesting to see, whether one can extract F2 in a 



closed form directly from the recursion rule (60) without solving for F^"^. Or one could try 
to prove that (68) solves the recursion rule ( |60[ ) by explicit computation or some inductive 
method. We have attempted this, but the problem is more complicated and we will proceed 



with assumption (67) 



We plug this formula for F2 into the the generating function (66), which after some 
algebra can be brought into the form 

Mt) = - ^ + ^aa" + f ) a/T:^. (69) 



13 



We have recovered the standard generating function for the distribution of M. Now, the 
discussion goes along the usual lines, the resolvent is 



1 



a;(A) = ^</.i(l/A) = dx^, 

A J -2a A - X 



which yields for the distribution 



(70) 



(71) 



This is the polynomial deformation to the Wigner semicircle distribution, with the radius 



2a. Looking back at (57), the variable in the case of the unsealed matrix is x/a/J and we 
need to replace g — )• j'^g ^i.e. the final formula is 



-{- 



+ 45/0^ + 2gx^ \/4aV-a;2. 



(72) 



Therefore starting from the recurrence relation (60) we have been able to recover the 



result of the distribution of eigenvalues of the random matrix ensemble with weight (53). In 



the next section, we will generalize this approach to different observables of the interacting 



theory. After using the explicit formula (45) for / we see, that the expression (72) reduces 



to the previous result in [10], where a polynomial deformation of the Wigner distribution 
was obtained alsoj£] 



4.2 Eigenvalue distribution of the matrix B and the joint MB distribution 

As in the free case, the theory now includes new observables involving the matrix B = KM. 
In the free case, this matrix followed the same distribution as the underlying matrix M. 
Now, the situation is going to be different, since the interaction involves only the matrix 
M. 

We will discuss the distribution of eigenvalues of B and the joint distribution for M 
and B. These two are going to be connected, since contractions of B with matrix M in the 
interaction vertex are going to turn even pure B correlators into mixes MB ones. Define 



' ' a 



2b 



1 



1 



Tr 



Tr 



M 

77 



2b 







Tr 









a times 



Tr 



M 

77 



(73) 



-V 







Tr 







a times 



Tr 



O,con 

M 

77 




^To see this better, we should car ry the explicit factors of / in the calculation. Such calculation would 
yield an e xtra factor of f"^+^"- in ( 68 1. From the definition of 4>i (t) we can see that /™ rescales t by v7 and 
from (581 we can see that J^" rescales g hy f^. 

^To do this, one has to introduce a parameter e in front of the kinetic term, do the expansion of / in 
powers of e and take e = 1 at the end of the calculation. 
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and 



2b 



w. 



m.b 



a=0 



a=0 



a=0 



a=0 



Using the same approach as before, it is now quite easy to write down the recurrence rule 
for 



p=0 a=0 



(76) 



where the first term comes from the contraction of B with a vertex. And again, this holds 
only for 6 > 1. It is also not too difficult to write down the recursion rules for W^s. Here, 
we can obtain two different recursions, considering contraction of the first M matrix or the 
first B matrix. These two are 



m—2 a 



6-1 a 



(77) 



p=0 0=0 
holding for m > 1, a > 1 and 

6-2 a 



p=0 a=0 



m—1 a 



wr^' = E E Giwr^'-'-' + 7 E E ^s^^'-a''''"' + 47w^r-t 

p=0 a=0 p=0 a=0 



(78) 



for 6 > 1, a > 1. Note that for 6 = in the first and m = in the second we recover the 
recursion rules for i^'s and G's respectively. Also in these expression W^'^ are considered 
as initial values, given by the recursion rules of the free case from the previous sections. 

Setting m = 1 in the first of the recursion we get 

6-1 a 



A^w'j, = 7^^i'^^-^-f EE^K 



^2fe/=f2(6-l-fc) 
a ) 



(79) 



k=Oa=0 



where we have used the fact, that Wa''^ = is nonzero only for even b. Using this in the 
the recursion rule for G (76) we find 



6-1 a 



Gf = iWy'-' + (1 - 7^) E E G'jGlt'-'^ 

p=0 a=0 



and 



6-1 



G2b = 7^^1,26-1 + (1 - 7^) J2 G2pG2(b-l-p)- 

p=0 

We again define the generating function </>2(s) = Yl s^^G2n, which yields 



(80) 



(81) 



Ms)-l = {l-7')s'r2is) + is 



1,26-lS 



26-1 



(l-7^)sV^(s)+7s 



.6=0 



,(82) 
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where we have used the fact that Wi^b vanishes for even b. We see that this is very different 
from the equation for <?ii(t), as expected due to the different role of M and B in the 
interaction. The solution is given by 



1 - y/l - 4s^{l - + -fsWi 
2(1-72)52 



Ms) = 7^, 2W • 



From the condition (p2 — ^ (po in the limit of g — )• we recover that in this limit Wi — )• 'yscpQ. 
We will prove that this is indeed the case shortly. 

In the following, we will drop the argument of 0i and (p2- It will be understood that the 
former is always functions of t and the latter function of s. 

From the recursion rules for T^'s, we can derive an equation for the generating function 
of the two point distribution 

^it,s) = J2t""s'W^,b- (84) 

m.b 



To do this, we express (77) as 

m-2 6-1 

p=0 p=0 
- 4.gWm+2,b- 

The terms with a subscript will cancel, since they follow the modified relations for the 
free quantities ([6]). Continuing the procedure we arrive at 

0(t, s) - 02 = t^</'i0(t, s) + ^ts^2Ht, s)-^ [Ht, s)-<P2- tWiis) - t'^W2is)] , (86) 

where we have denoted W2{s) = Ylb^o^'^fiS^ ■ ^ similar equation can be derived from the 
second recurrence rule 

</.(*, s) - = s^(l)2Ht, s) + ^tscPi^Pit, s)-^ [</.(t, s)-cl)2- tWiis) - t^W2{s)] . (87) 

Before, we have used multiple equations for the generating function as a consistency check. 
However by now we trust our procedure enough to using these two together to reduce the 
number of unknown functions to one. Doing this, we get the final formula for (f)(t, s) 



t,s) 



(t>2 {t(t)l - IScf)--, 



t — 'ys(j)2 + jstWi 

So to compute 02 (s) and 0(t, s) we need to specify the function Wi{s). 

Let us stress here, that the coefficient 7 in these expressions is the very same correlation 
as in expressions of section [2| This shows, that even though we consider different matrix 
ensemble, it still knows about the underlying fuzzy sphere it wast built on, which is therefore 
encoded solely in the kinetic term of ([T]). 
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Figure 2: Triangles represent matrices B. With only one MB contraction, the diagram 
splits into two parts of matrices B and contractions among M's. 



4.3 Leading order Wi{s) in 7 

In the case of F^™', we have been able to guess the solution of the recursion rule. However for 
the case of Wa''^^~^ the situation is more involved and no easy guess is possible. Moreover, 
simple analysis of the integer factorization of the first couple of terms shows, that such 
guess might be very difficult and the explicit formula much more complicated than a simple 
product of factorials. 

Therefore we will compute only in the leading order in 7 to get the first non- 

trivial contribution due to the MB contraction. In such case, the contributing diagrams 
have only one AdB contraction, leading to one factor of 7. Let us stress that the results 
we will obtain are still exact in g and we make no assumption about the magnitude of the 
coupling. 

As figure 4.3 illustrates, we get the following expression for the diagrams with one MB 
contraction 



p=0 



^^1,26+1 = ^ ^ CpCb-p = 7Cb+i = 7F0 Cb+1. This then yields 

00 

W,^2b+i = Y.(-9rwy'^' = 7C6+ii^2 (90) 



where we have used the recursion rule (60). This holds for a >. For a = we simply have 



a=0 



and 



W, = ^.^1,, = Y.s''+'W,,2b+^ = ^F,^^^^^^ = ^s4>l{s)F,. (91) 

b=0 b=0 
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We use this result in the equation (83 ) to obtain the leading contribution to the distribution 

02 



MS 



—— = 0o(s) +7 (^2 - ^m{s)- ^-27^ ■ 



(92) 



In the limit of g' — )• 0, F2 becomes 1, the whole extra contribution vanishes and 02(5) 



Expanding (88) up to second order in 7 yields 

^i(*) 



(t){t,s) = 0o(s)'Ai(*) +7'S</'o(s)- 



1 



+ 



+ Y 



Mt){F2-i)Ms) 



t 

S^lis) {l-F2t^)M t) - 1 2 

1 - smis) i2 



(93) 



We now need to invert these to obtain the distributions. We can use the standard 
approach of defining the resolvent ( 70 ) and then using the discontinuity equatiorj^ 

1 



p{x) 



2-Ki 



[ljj{x + ie) — u}{x — ie)] 



(94) 



But it is easier to use a fact mentioned in the appendix |A] that the non-singular part of 
(f){t)/t'^ generates the distribution x^p{x), where p{x) is generated by (pit). This way we 
find out that 



(95) 



generates y^(l — y^)po{y) and similarly for ((1 — F2t^)4>i{t) — l)/t^ ■ After some algebra, we 
find final result for the distribution of the eigenvalues of matrix B and the joint distribution 
of M and B, valid up to second order in the correlation 7 to be 

,2 



P2{y) 
p{x,y) 



Po{y) l + 7'(l-i^2) 



y 



y^ 



Pi{x)p2{y)(l + ixy + j X {I-F2X )y {l-y 



(96) 



As in the case of the distribution of M, we have to change g — gf'^,x — >• x/\/J,y — )• y^/g 
to get the distributions of the unsealed matrices. Note that both of these become the 
appropriate free expressions but there is something new in the second formula. One could 
guess that the interacting result for p{x, y) would be just the free case with the one matrix 
marginals replaced by the interacting expressions. The extra factor of F2 shows that this is 
not the case. 



5 Conclusions and Outlook 



The main results we present are twofold. First, we have shown that it is possible to compute 
distributions of more general observables in the matrix ensemble discussed in [8j . Also, we 

^For more details, see e.g. [13| 
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show how rescahng of the mass puts ah the terms on the same footing in the large A'^ 
hmit. Second, we have shown that the matrix ensemble corresponding to the interacting 
scalar field on the fuzzy sphere does reproduce the distribution of eigenvalues of the model 
without the kinetic term, which is the polynomially deformed Wigner semicircle. The 
radius of the distribution gets renormalized by a factor of as was the case for the free 
field ensemble and this therefore seems to be a generic feature. We have obtained the 
equations for the distribution of matrix B and the joint distribution of M and B and we 
have solved these up to an unknown function W\. We have then computed the first non- 
trivial contribution, which enabled us to compute these distributions to the second order in 
7. The results we have obtained show, that the interaction introduces a novel features into 
the joint distribution of M and B beyond changing the one matrix marginals. 

There are two points that deserve further treatment. First of all, the proof that the 



formula ( |68| does indeed solve the recursion rule (60), or equivalently recovering the ex- 
pression ( |67[ ) for F2 directly from the recursion. Also more systematic treatment of the 7 
expansion of ^2('S) and s) is needed and could lead to better understanding and possibly 
complete solution for VFi(s). 

With the results presented and possibly some future progress there are several issues 
to explore. The presented results treat the kinetic term non-perturbatively, therefore it 
is going to be interesting to analyze the phase diagram of the theory and to compare the 
findings with previous numerical work |14j . It has been also shown that a modification of 
the kinetic term can remove the tadpole diagrams responsible for the UV/IR mixing |15| . 
Our method is well suited for such modification and one will be able to study signatures of 
this removal in phase diagram of the theory. 

Acknow^ledgments. I am grateful to V.P. Nair, A. P. Polychronakos and O. Budac for 
a lot of fruitful discussions. This work was supported by U.S. National Science Foundation 
grant PHY-0855515 and by a PSC-CUNY grant. 



A Proof of the formula (29) 



We will prove the one dimensional case. The generalization to more dimensions is then 
straightforward. The presented proof is a modification of approach used in [8]. 

Let us have a generating function of the form 

^{t) = f{t^{tmt), (97) 

where (l>{t) = 2/(1 + \/l — t^) is the generating function for the Wigner semicircle distribu- 
tion. We expand this as 

00 

$(f) = ^a„t>"+i(t). (98) 

n=0 
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Using the explicit formula for <j){t) we obtain 



A:=0 ^ ' / V / fc=o 



n 



2/c + l/ 



(l-t^)^ (99) 



Now follows a crucial observation. If a generating function ^(t) corresponds to the dis- 
tribution p{x), then the non-singular part of g{t)/t"' generates x^p{x) for any n > 0. In 
previous expression, the second term contains only singular terms. However the left hand 
side is clearly non-singular and therefore the second term cancels the singular part of the 
first term but does not contribute otherwise. After resuming the first term of (99), we see, 
that t'^(j)^^^{t) generates the following distribution 



1 



2./l-^ 



1 



n+l 



1 



71+1 



2" 

ITT 



+ i\/ 1 — 



n+l 



X — I 



a/I — 



n+l 



Since |x| < 1, we can write x = cosO for some 6 G [0,7r] and x it i^/Y 



yields 

Pn{x) = 

And the final distribution is then 



2n ^i{n+l)9 _ 2n^-i{n+l)B 

in 



(100) 
5±*^. This 

(101) 



n=0 



X} = — I e*V(2e*'^) - e-*V(2e 



(102) 



which is the desired formula (30) for the one dimensional case. 



B Explicit formulas for three and four matrix distributions 

The factor p3{x,y, z) multiplying p{x)p{y)p{z) in formula (31) is given by fraction with the 
following numerator 

1 - 9l2 - 5l3 + 512513 - 912913923 + 9l29l3923 + 5'12S'l35'23 - 9l29l3923 " 523 + 5l2523 + 
+513523 - 5l25l3523 + 512513523 " 5l25l3523 " 5l25l3523 + 5l25l3523 " '^912913^'^ + 
+45125135232;^ + 45125135232;^ " ^912913923^^ + ^9l29l3Xy " ^^913923X9 - 45125135232:2/ + 
+^9l29l3923Xy + 45i25i32;y - 4:912913923^9 - 45125135232^5 + 45125135232^5 + 45125135235^ - 

45125135235^ - 45?25235^ + 45i25i35235^ + 49u9i3xz - Agi2g23xz - Agi2gl^g23xz 
+4gl29i3923xz + 45135I3XZ - 45125135232:2 - 4gl29i3923xz + Agi2gl^g%xz - 4gi2gi3yz 

+451252352; + 45135232/2 - 45^251352352 - 451251352352 + 451251352352 
+451251352352 - 45i25l352352 + 45125135232^ - 45^25135232^ - 45135232:^ + 45i25?35232^ 
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and denominator 

o222|^422 r)42|^242 o442,4 o24,44 o24 
-87i27l3723 + 47l27l3723 - 2713723 + 47l27l3723 - 27l27l3723 + 723 - 2712723 + 712723 - 2713723 

+47?27?3 713-27^2 7^3 71^+7^3 7i3-27f2 7^3 7^3+7^2 7^3 7i3+47?2a;^+47L3-'c^-167?27?3-^^+47^27?3a^^ 
+47i27i3a;^-87i27i3a:;^-87?3723a:^^+327L7i3723a:^^-87i27i3723a^^-87i27i3723aJ^+47i^ 

+47137233^^ - 167127137232?^ +47l27l3723^^ +47127137233^^ + 167^27133^^- 327f27l3723 

-4ji2xy - Ajf^xy + 87i27i3a;y + 87?27i3a?y - ^'H2ltiXy - ^Julfs^V + ^-fulh^V + 87?2723a?2/ 
-167i27?37i3a?y-167?27f3723a^?/+87i27i3723a^y+87?27^37l3a?2/-47i27l3a^y-47?27l3a?y+87i27?37t3^2/ 

+87i27?3723^'y - 47127^3 7|3a;y - 47^27^3713x7/ - 16712713^'^y " 167^2 7i33?^y + 327i27i37232;^y 
+327?27i3723a;^y-167i27f3723a;^y-167?27f3723a;^y+47i2y^-87i27i3y^+47i27i3y^+47l3y^-167i2723y^ 
+47^2 7l3y^-87?37l3y^+327f27f37i3y^-87f27?37l3y^+47f37l32/^-167f27f3723y^ +47^27^3 7^3^^ 

, A 2 4 2 o2 2 4 2,^2 4 4 2iic2 2 22iic2 2 22iit;2 2 22 nc!2 2 2 22, 

+4712723^ -87x27x3723^ +47i27i3723y +167127132? y +16712723^? y +16713723^^ y -96712713723^; y + 

ic4 2 2 22iic2 4 2 22iic2 2 4 22ic,i2 2 2 42 ic 2 3 ia 3 2 3 

16712713723^? y +167127137232? y +16712713723^? y +64712713723^; y -167127232^2/ -167i27232;y 

+327i27?37l3a?y^+327?27f37l32?y^-167i27t37i32?y^-167?27t37i32;y^-647i27f37i32;^y^-647?27i^ 

+167i2723/-327?27i3723/+167i27i3723y^+647?27i37232;^/-47i3a?2;+87i27i3a?2;-47f27i3a?2;-47?3a;z 

+87i27i32;z-47^27i32;2;+87i37232;2-167?27i37232?2+87i27i3723a?2+87?3723a?^-167i27i37232;2 

+87127137232?^; - 47137232^2: + 87?27l37232;2: " 47^271 3 7132^2; - 47?37232;2 + 87?27X37232;2 

-47f27?37232;2-167?27x32;^2-167i227?32;^2+327f27x37i32;^2+327f27?37i32;"'^2-167i227i37l32;^2 

-167l27l37232;^2-4723y2+87i2723y2-47i2723y2+87i3723y2-167i27l3723y2+87l27l37232/2;-47i3723j/2 
+87127^3723^2 - 47^27^3723^2; - 47132/2; + 87^2723^2 - 47^27232/2 + SlfslisV^ - 16712713723^2 

+87i27i3 713^2 - "^itslhyz + 87?27i3723y2 - 47^2713723^2 + 167i27i32;^2/2; + 167?27i32;^y2 
+167i27i32?^y2+167i27i32;^y2;-167i2723a;^y2-167f3723a?^y2;+647i27i3723a;^y2;-167i27i3723a:;^y2; 
- 167f27i3723a?^y2;- 327i27i37232?^y2;- 327^27i37232?^y2; - 32712713713X^^2; - 327^27i37232?^y2 

-167l27232?^2/2;-167f37f3X^J/2;+647^27l37232?^2/2-167t27l37232?^y2;-167f27l37232?^y2;+167i27l37232;^y2; 

+167i27i37232;^y2+167i27?37232;^y2+167?27i37232;^y2-647^27i3723a;^y2-647f27i37232;'^y2 
-167f27i3a;y^2;-167i27i3xy^z+167i2723a;2/^2+167f2723a;y^2;-327i27i3723a:y^2-327f27i3723a?y^2 
+167i27i37232;y^2+167i2 7^3 7232;y^2-167i37|3xy2z+647i27i37232;y^2-167i27i37232;y^2-167i3 713x2/^2; 
+647i27i37232;y^2-167f27i37232;y^2+167i27232;y^2+167f27i32?y^2-327i27f37i3xy^z-327i'27i37232;y^2; 

+ 167l27l37232;y^2+167i27]'37232;y^2-167i27l37232;y^2-167i27i37|3Xy^Z+647i27i37232;^y^Z 

+647?27?3723a?^y^2-647i227i37l32?^y^2;-647i227?37l32?^y^2;+647i27?37i32?^y^2;+647i327i37i3;^^y^2; 
-167i2723y^2;+327f27f3723y^2;-167?27i3723y^2;- 167^27232/^ 2;+327f27f3723y^2;- 167^27137232/^2; 
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-647f27f37232;^y^z+647i27i37l32;^y^^+647?27i37l32;^y^^+647i27?37l32;^y^^+647?27?37i32;^y^^ 
-647^2 7?37i3^^y^^ - 64712 7i37i32;y'^^ - ^^lulislh^v"^ ^ + 47^3^^ " 871^2 7?3^^ + '^lulh^'^ 

+4723^^ - 8712723^^ + ^l\2llz^^ - 167?37l3^^ + 327?27l37l3^^ - 167127137232^^ + 4713723^^ 

-87?27l37l3^^ + 47^27^3713^^ + 47?3723^^ " 87?27?37l3^^ + ^Il2ll?.l2-i^'^ + 167^27?3^^^^ 

X Mi 2 2 22iic2 2 22 nc 2 2 2 22|ic4 2 2 22,1^2 4 2 22,1^2 2 4 22 
+ 16712723^; Z +16713723^^ ^ -96712713723^^ ^ +16712713723^^ ^ +16712713723^^ ^ +167127137232^ ^ 

+647^2 7i37232;^^^-167i27i3a^y2:^-167i27i3a;y2;^+167i37232;y^;^-327f27i3723a::yz^+167f27i3723a:y^^ 

+167^3 723xyz^-327i27i3723a:^y2;^+167i27i3723a;y^^-167i27232;y2:^-167i2723a^y^^+647i27i3723^y^^ 

+64712 7i3723a;y2^-167i27i3723a;y^^-167i27i3723a^y^^+167i3723a;y^^-327i27i37232;y^^+167i27i^ 

+167?3723a;y^^ - 327i^27?37f3a^y2:^ + 167i27i3723a^y^^ - 16712713723^^^^^ - 16712713723^^^^^ 

+647f27i37232;^yz^+647f27?37232;^yz^-647i27?37i32;^y^^-647?27?37i32;^y^^+647?27i37i32;^y^^+ 

647i'27?37l3^'y^' + 167?27?3y'^' + 167?27i3y'^' + 167?37l3y'^' " 967?27?37l3y'^'+ 

167i27i3723y^^^+167i27i3723J'^^^+167i27i3723y^2;2-647i27i3723a;^y^2^-647j'27i3723a;^y^^^- 
647i27?3723x2y^2^-647i327?3723a;^y^2^+1287i227f37|3x2y2^2_g4^^^^^^^3_^^2y2^2_g4^3^ 

647i27?37i3^S^^^-647?27?37i3a;^y^^^+647f27i37232;y^z^+647f27?37232;y^z^-647i27?37l32;2/^^^- 
647?27?37l3^y^^^ + 647i227i37|3xy3^2 ^ U-i\^-il.i-il.^xy^ + 6471227123723^^^2 _ x^^^^f^^xz"^ 

+327?27i3723a;2^-167f27i37l3a;^^-167i3723a^^^+327?27?3723a;^^-167i27i37l3a;^'^-647?27i3723^^^^- 
647?27i37232;^^^-167i3723y^;^+327i22 7^3 723y2^-167f27i3723y2;^-167?37l3yz^+327i22 7137^3^2^- 

167i27i3723y^^-647i27i37232;^y2;^+647i27i37232;^y^^+647i27i37232;^y^^+647i27i37232;^y^^+ 

647?27?3723^^y^^-647?27?37l3^^y^^+647i27?3723xy2z3+647?27?3723a;y^^^-647i227i37l3a^y^^^- 

647f27?37l3^y^^^+647i27?37l3^y^^^+647?27?37l3^y^^^-647?27i3723y^^^-647i227?37l3y^^^+167?37l3^^ 
-327?27i3723^^+167i27i37l3^^+647f27i3723^^^^-647i27i37233;y^'^-647?27i37232;2/2'^+647227i237|3y2^4 

The factor p4{x,y, z,w) multiplying p{x) p{y) p{z) p{w) in formula pS] ) is given by 

(1 - 72)3 [1 + ^6 _ ^2(^]^ _^ ^2-^(j_ _^ 4^,^ _^ 4y^^) _^ 4^1^^ _^ _^ 

denominator 

with the denominator 

(l+72)^-47(l+7^^)(x+z)(i(;+y)872(l+7^2^ x^+y^+z^+w^+2{xy+zw){xw+yz)+Axyzw 

-47^(1 + 7^°)(x + z){w + y) - 5 + 4(^2 + ^2 + ^2 + w'^) + A{xz + wy) + Idxyzw 

+167^(1+7^) -Six"^ +y'^ +z'^+w'^) + {x'^+y'^+z'^+w'^)+3ix'^ +z'^){y'^ +w'^)+4{x'^z'^ +y'^w'^) 
—Axyzw — 2{xy + zw){xw + yz) + 8xyzw{x'^ + y'^ + z"^ + w"^) + 4(^2 + z^){y'^ + w'^){xz + wy) 

+\Qxyzw{xz + wy) + KSx'^y^ z^w'^ 
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87^(1 + 7^) 



47^(1 + 7 + z){w + y) 9 - 12(j;^ + + + + xz + wy) + 16(x^ + z^){y^ + w^) 
+16(x^2;^ + y'^w'^) + 16x2;(x^ + z^) + 16yi(;(y^ + 'uP') + 48xyzw + 6Axyzw{xz + wy) 
15(x^ + + + w^) - 8(x^ + / + z'' + ti;"^) - 2{xy + + yz) 

-24(x^ + 2;^)(y^ + u;^) - 32(x^z^ + y^iD^) - 4xyzw + 32(2;^/2;^ + perm) + 16x^z^(x^ + z^) 
+16/t(;2(y2 + w^) + 8(j;2 + z'^){y^ + u;^) + 8(x'' + z^)(y2 + w^) 
+32xy2;i(;(a;z(x^ + z'^) + yw{y'^ + w'^)) + 2>2xyzw{x^ + z^){y^ + 
+64:xyzw{x'^ z"^ + y'^up') + 6Axyzw{xy + zw){xw + yz) + 32(x^y^z^ + cycl) + 12d>x'^ y^ z^nP 

-Wil + 7^)(x + z){w + y) - 5 + 8(x^ + y^ + + u;^) + S{xz + yw) 
— 16(2;^ + z^){y'^ + nP) — \Qxz{x^ + z^ + xz) — 16yw(y^ + w'^ + yw) — 64xyzw 
+64:xyzw{x'^ + y^ + z^ + uP) + 128xyzw{xz + yw) + 6A{x^z'^ + y'^w^) + 64(x^y^z^ + perm) 

327^ - 5{x'^ + y^ + + w"^) + 3(x^ + y^ + + w^) + 9(2;2 + z'^){y'^ + ti;^) 

+2(xy + zw){xw + yz) + 12(x^z^ + y^w^) + Axyzw - 8x'^z'^{x^ + z^) - 8y^tt;^(y^ + z^) 
— 4(x^ + z'^){y'^ + w^) — 4(x'^ + z'^)(y^ + w^) — 8xyzw{x'^ + y^ + z^ + w^) — 16xyzii;(xz + yw) 

-16(x^y^z^ + perm) - 4(x^ + z^)(y^ + w^){xz + -wy) + 8(x'^z^ + yW) + 64x^y^z^zi;^ 
+32xyzty(x^z^ + y'^w'^) + \Qxyzw{x^ + z'^){y'^ + w"^) + lQxyzw{xz{x^ + z^) + yw{y'^ + w"^)) 
+32xyzw{xy + zu;)(xtt; + yz) + 16(x^y^z^ + cycl) + 8(x^y^z^ + x'^y^z'^ + x^y^z^ + perm)) 

Note that all the polynomials have the desired symmetry x ■h^ z,y ■h^ w, (x, y) o (z, tf). 
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